Relativistic quantum chemistry on quantum computers 
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Last years witnessed a remarkable interest in application of quantum computing for solving prob- 
lems in quantum chemistry more efhciently than classical computers allow. Very recently, even 
first proof-of-principle experimental realizations have been reported. However, so far only the non- 
relativistic regime (i.e. Schroedinger equation) has been explored, while it is well known that 
relativistic effects can be very important in chemistry. In this communication we present the first 
quantum algorithm for relativistic computations of molecular energies. We show how to efficiently 
solve the eigenproblem of the Dirac-Coulomb Hamiltonian on a quantum computer and demonstrate 
the functionality of the proposed procedure by numerical simulations of computations of the spin- 
orbit splitting in the SbH molecule. Finally, we propose quantum circuits with 3 qubits and 9 or 
10 CNOTs, which implement a proof-of-principle relativistic quantum chemical calculation for this 
molecule and might be suitable for an experimental realization. 



Quantum computing [l| is one of the fastest growing 
fields of computer science nowadays. Recent huge in- 
terest in this interdisciplinary field has been fostered by 
the prospects of solving certain types of problems more 
effectively than in the classical setting [3, [l] . The promi- 
nent example is the integer factorization problem where 
quantum computing offers an exponential speedup over 
its classical counterpart 0] . But it is not only cryptogra- 
phy that can benefit from quantum computers. As was 
first proposed by R. Fcynman quantum computers 
could in principle be used for efficient simulation of an- 
other quantum system. This idea, which employs map- 
ping of the Hilbert space of a studied system onto the 
Hilbert space of a register of quantum bits (qubits), both 
of them being exponentially large, can in fact be adopted 
also in quantum chemistry. 

Several papers using this idea and dealing with the in- 
terconnection of quantum chemistry and quantum com- 
puting have appeared in recent years. These cover: cal- 
culations of thermal rate constants of chemical reactions 
0, non- relativistic energy calculations quantum 
chemical dynamics [l^l, calculations of molecular prop- 
erties [ill , initial state preparation [l2|, |l3| , and also first 



proof-of-principle experimental realizations 14|-|17|. An 



interested reader can find a comprehensive review in [l8l | . 

An efficient (polynoniially scaling) algorithm for cal- 
culations of non-relativistic molecular energies, that em- 
ploys the phase estimation algorithm (PEA) of Abrams 
and Lloyd [l^ . was proposed in the pioneering work by 
Aspuru-Guzik, et al. 0] . When the ideas of measurement 
based quantum computing are adopted the phase 
estimation algorithm can be formulated in an iterative 
manner [iterative phase estimation (IPEA)] with only one 
read-out qubit If the phase </> (0 < (/> < 1), which is 

directly related to the desired energy is expressed in 



the binary form: (j) ~ 0-4'i4'2 ■ ■ 4'i — {0, 1}, one bit of (j) 
is measured on the read-out qubit at each iteration step. 
The algorithm is iterated backwards from the least signif- 
icant bits of (f> to the most significant ones, where the A:-th 
iteration is shown in Figure [TJ Not to confuse the reader, 
H in the exponential denotes the Hamiltonian operator, 
whereas H (in a box) denotes the standard single-qubit 
Hadamard gate. |?/'systGm) represents the part of a quan- 
tum register that encodes the wave function of a studied 
system, is a z-rotation gate whose angle depends 
on the results of the previously measured bits [Sl, |9[ , and 
parameter r ensures that < (j) < 1. The PEA always 
needs an initial guess of the wave function correspond- 
ing to the desired energy. This can be either the result of 
some approximate, polynomially scaling ab initio method 
or as originally proposed by Aspuru-Guzik, et al. 
the exact state or its approximation prepared by the 
adiabatic state preparation (ASP) method. 
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Figure 1. The fc-th iteration of the iterative phase estimation 
algorithm (IPEA). The feedback angle ujk depends on the 
previously measured bits. 



It is a well known fact that an accurate description of 
molecules with heavy elements requires adequate treat- 
ment of relativistic effects [2lj. The most rigorous ap- 
proach [besides the quantum electrodynamics (QED) 
which is presently not feasible for quantum chemical pur- 
poses] is the four component (4c) no-pair formalism. Our 
work is based on the 4c electronic Dirac-Coulomb Hamil- 
tonian (DCH) in the form 
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Dirac matrices appearing in the one-electron part are de- 
fined as 
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the former in terms of the Pauh spin matrices <t. The 
DCH is known to cover the major part of the spin-orbit 
interaction and also scalar relativistic effects. Using this 
type of Hamiltonian represents no loss of generality for 
our purposes, since a transition to the Dirac-Coulomb- 
Breit Hamiltonian 221 and the inclusion of the corre- 



sponding integrals requires a classically polynomial ef- 
fort. 

We also adopt the no-pair approximation (NPA), 
widely used in relativistic quantum chemistry [22| . in 
which the N-particlc basis of Slater determinants is con- 
structed from positive-energy bispinors only. For a more 
detailed discussion about the DCH and approximations 
employed in the relativistic quantum chemistry, see the 
Supplementary Information. 

The use of a 4c relativistic formalism brings in three 
major computational difficulties compared to the non- 
relativistic case: (1) working with 4c orbitals (bispinors), 
(2) complex algebra when molecular symmetry is low, 
and (3) rather large Hamiltonian matrix eigenvalue prob- 
lems [due to larger mixing of states than in the non- 
relativistic (NR) case] . The central objective of this work 
is to address these problems in regard of an application 
of a quantum computer and the extension of the quan- 
tum full configuration interaction (qFCI) method to the 
relativistic regime. 

We will start the description of the algorithm with a 
mapping of the relativistic quantum chemical wave func- 
tion onto a quantum register. The simplest (scalable) 
NR approach, the direct mapping (DM) , assigns each 
spin orbital one qubit (|0) = unoccupied, |1) = occu- 
pied). The relativistic case is similar due to the NPA. 
Moreover, because of the time-reversal symmetry of the 
Dirac equation, bispinors occur in degenerate Kramers 
pairs [22| denoted A and B (analogy to a and (3 spin in 
NR treatment) and the relativistic DM thus looks like: 
one qubit for bispinor A and one for B. The 4c charac- 
ter of molecular bispinors therefore does not complicate 
the approach substantially [note that as in the NR case, 
the Hartree-Fock (HF) calculation is done on a classical 
computer and only the exponentially scaling FCI on a 
quantum one] . 

The DM is known to be not optimal as it maps the 
whole Fock space of the system on the Hilbert space of 



qubits. For this reason, compact mappings from a sub- 
space of fixed-electron-number and spin- or symmetry- 
adapted wave functions have been proposed 0, . How- 
ever, general factorization schemes [i.e. algorithms to 
systematically generate quantum circuit implementing 
exp(iT/f)] for these mappings have not been discovered 
yet. In the relativistic case, the most convenient compact 
mapping is based on a subspace of symmetry-adapted 
functions employing the double group symmetry. 

Assuming the NPA and the empty Dirac picture, the 
relativistic Hamiltonian has the same second quantized 
structure as the NR one 
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Here hpq and 



■s denote one- and two-electron integrals 
that are in contrast to NR ones in general complex. This 
is in fact no difficulty for a quantum computer, since our 
working environment is a complex vector space of qubits 
anyway and we do the exponential of a complex matrix 
even if the Hamiltonian is real (see Figure [T|) . After the 
decomposition of the unitary propagator [exp(iriJ)] to 
elementary quantum gates_(in case of DM) using the 
Jordan- Wigner transform [23|, one can see that complex 
molecular integrals require twice as many gates compared 
to real ones [8|, while complex arithmetic on a classical 
computer requires four times more operations. 

The last of the aforementioned difficulties of the 4c 
formalism is the size of a Hamiltonian matrix eigenvalue 
problem. This can be inferred from the observation that a 
significant larger number of integrals in the Hamiltonian 
^ will be non-zero due to the lowering of symmetry in- 
duced by spin-orbit interaction. The loss of spin symme- 
try can to some extent be alleviated by consideration of 
time reversal symmetry. In the Kramers-restricted (KR) 
approach employed in this work the second-quantized 
Hamiltonian ([3]) is expressed in terms of a basis of 
Kramers pairs, that is, orbital pairs (j) and 4> connected 
by time reversal. Determinants may be characterized by 
a pseudo-quantum number Mk = 1/2{Na — Nb), re- 
flecting the different number of unbarred A''^^ and barred 
Nb bispinors. In the non-relativistic limit the Kramers 
pairs can be aligned with spin partners such that Mis- 
becomes identical to Ms- However, contrary to the NR 
limit, determinants with different Mk can mix in the 
presence of spin-orbit interaction. It can be shown (see 
Supplementary Information) that the ratio between di- 
mensions of relativistic and non-relativistic Hamiltonian 
matrices scales as 0(m^/^) in the number of molecular 
orbitals/bispinors. 

When employing the DM on a quantum computer, this 
problem does not occur, since the Hamiltonian ([3]) then 
implicitly works with all possible values of AIk. The 
scaling of the relativistic qFCI method is therefore the 
same as the NR one, namely ©(m^) d, 14| , where m is 
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GAS 


Min. el. 


Max. el. 


Shell types 


I 





4 


0"l/2, 1"i/2 


II 


2 


4 


1"3/2 


III 


4 


4 


(jjyj, 43 virtual Kramers pairs 



Table I. GAS and occupation constraints for SbH X and 
A 1 states CI calculations. The minimum and maximum num- 
ber of electrons are accumulated values - apply to this and all 
preceding GA spaces. 

the number of molecular orbitals (bispinors). 

For numerical tests of the algorithm, we have cho- 
sen the SbH molecule whose non-relativistic ground state 
•^E" splits due to spin-orbit effects into X 0+ and A 1. 
In the approximate A^^ notation, these states arc dom- 
inated by '^'^f2'^'i/2^^/2 ^'^'^ '''1/2 ""1/2 ""3/2 configurations. 
The splitting is truly of "molecular nature" as it disap- 
pears for dissociated atoms. Its experimental value is 
A£;so = 654.97 cmT^ (H. 

In all our simulations, we used the Dyall triplc-zeta 
-I- valence correlating functions, total 28s 21p 15d If for 
Sb and cc-pVTZ (from EMSL basis set library) for H. 
We, of course, could not manage to simulate the FCI cal- 
culations with all electrons in such a large basis. We in- 
stead simulated general active space (GAS) KRCI com- 
putations [25| with the occupation constraints shown in 
Table U giving rise to CI spaces of approximately 29500 
determinants. For a balanced description of both states, 
we optimized the spinors taking an average energy ex- 
pression (2 electrons in 2 Kramers pairs 1^1/2, ''■3/2 )■ 
We worked solely with a compact mapping employing 
the double-group symmetry (C|.jj) and exponential of a 
Hamiltonian was simulated as an n-qubit gate (similarly 
as in d, 0,0). We used the DIRAC program for cal- 
culations of Hamiltonian matrices. The nuclear potential 
4>nuc was generated by finite nuclei using Gaussian charge 
distributions with exponents chosen according to Ref. |27 . 
Simulations of qFCI computations were performed with 
our own C++ code 0. We ran 17 iterations of the IPEA 
with the difference between max. and min. expected en- 
ergies equal to 0.5 -E^ We also did not count the least sig- 
nificant binary digit of the phase to the total success 
probability (for more details of the algorithm, we refer 
the reader to our preceding paper Q). This procedure 
corresponds to the final energy precision w3.81 x 10^^ 
Eh. 

Simulated potential energy curves of both states are 
shown in Figure O Based on our KRCI setup we obtain a 
vertical AEso of 617 cm~^. Success probabilities (SPs) 
of the algorithm with HF initial guesses (o'i/2'''i/2'''3/2 

for A 1 one) are 



for the X 0+ state and f^^2'''i/2''''3/2 
presented in Figure [S] They correspond to the IPEA 
with the second part of a quantum register (encoding the 
relativistic quantum chemical wave function) maintained 
during all iterations (in [§] denoted as version A). In 
this case, SPs always lie in the interval |(V'init|V'exact)P • 
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Figure 2. (Color online) Simulated potential energy curves 
of ground (0"*") and excited (1) states of SbH, and spin-orbit 
energy splitting. Absolute energies are shifted by 6481 E^- 
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Figure 3. (Color online) SbH ground (0"'") and excited (1) 
state qFCI success probabilities (SPs) corresponding to HF 
initial guesses. 

(0.81, 1) 0. Ground state SPs confirm that relativistic 
states have, due to near degeneracies caused by the spin- 
orbit coupling, often a stronger multireference character 
than non-relativistic ones. The upper bound of the SP 
is less than 0.7 even for the equilibrium geometry and 
HF initial guesses can in fact be safely used (SP > 0.5, 
amplification of SP by repetitions) only up to 4.8 oq. The 
SPs of the A I state are higher and HF initial guesses can 
be in a noise- free environment used up to 6 oq. 

The difficulty connected with a low success probability 
for the X 0+ state at longer distances can be overcome by 
the ASP method 0. In this approach, one slowly varies 
the Hamiltonian of a quantum register, starting with a 
trivial one with a known eigenstate and ending with the 
final exact one in a following simple way 



H=il- s)Hi, 







1. 



(4) 



If the change is slow enough (depending on the gap be- 
tween the ground and the first excited state), the register 
remains in its ground state according to the adiabatic 
theorem (28j . In our relativistic example, analogously 
to the non-relativistic one m , iJinit is defined to have all 
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matrix elements equal to zero, except Hn, which is equal 
to the (Dirac-)HF energy. 
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Figure 4. (Color online) Adiabatic state preparation (ASP) 
of the SbH ground state (0^) for different internuclear dis- 
tances. Solid lines correspond to qFCI success probabili- 
ties, KV'ASpl'/'oxact)!^ ■ (0.81, 1) interval is colored. 1000 hE^'^ 
« 10-" s. 

We simulated X 0+ qFCI computations with adiabat- 
ically prepared states for different internuclear distances; 
results are shown in Figured In this case, for computa- 
tional reasons, we employed complete active space (CAS) 
KRCI method with a CAS composed of 2 electrons in 
the highest occupied (7^1/2) and 45 lowest unoccupied 
Kramers pairs (corresponds to 2116 determinants). It 
can be seen that for t = 1000 fiE^^ , the upper bound of 
the SP goes safely to unity even for r = 8 ao- 

Recently, there appeared two papers presenting the 
first physical implementations of non-relativistic qFCI 

15| quantum 



IJ and NMR 



computations on optical 
computers. Correspondingly, we would like to propose 
two candidates for the first relativistic computations on 
real quantum computers. Our proposals represent "dig- 
ital (circuit-based) quantum simulations" (DQS) as de- 
fined by Bulata and Nori [2^. Conceptually different 
are "analogue quantum simulations" (AQS), where the 
evolution of a studied quantum system is mapped to be 
simulated onto the controlled evolution of the quantum 
simulator. Recently, Gerritsma et al. used this approach 
for the proof-of-principle simulation of a one-dimensional 
Dirac equation with a single trapped ion (soj . 

Both of our examples represent calculations of SbH 
■^I]- ground state spin-orbit splitting. Since one has to 
employ rather large basis sets (triple-^ quality) to get a 
meaningful result, they again are not true FCI calcula- 
tions, but FCI calculations in a limited CAS. The first 
one corresponds to a CAS composed of 2 electrons in 
the highest occupied {Tri/2) and the lowest unoccupied 
(773/2) Kramers pairs [CAS(2,2)]. After the factorization 
of a Hamiltonian according to the quantum number 
and taking into account only one of the two degenerate 
2;-projections of f2 (for = 1), the size of the CI space is 
2 for the ground state (0"^) and 1 for the excited state (1). 
The excited state is therefore trivial and the calculation 



of the ground state is in fact a complete analogue of the 
already mentioned NR computations [l4,[l5|, because it 
needs just one qubit for the wave function (2 in total). 
The controlled singlc-qubit gate can be decomposed us- 
ing 2 controlled NOTs (CNOTs) Calculations with 
this active space yield an AEgo = 509 cm"^ computed 
at the experimental equilibrium bond distance of 3.255 
ao- 

The second example represents a 3-qubit experiment (2 
qubits for the wave function) and employs a CAS com- 
posed of 4 electrons in the cr\i2'^\l2''^'il2 Kramers pairs 
[CAS(4,3)]. It gives a better value of A£;so(518 cm"!) 
than CAS(2,3). After factorization, the CI space of the 
excited state has a dimension 3 and that of the ground 
state 5. Fortunately, near the equilibrium bond distance, 
the Hamiltonian matrix of the ground state is to a very 
good approximation block diagonal (ground state energy 
difference of the order /iSh), coupling only 3 configu- 
rations K/271-i/27r°/2, ^ll2Al2^ll2^ ^ud o\ / ^1^1 1 ^1^% i ^) . 

If we take into account only these configurations, both 
states can be encoded by two qubits. 

We used the Quantum Shannon Decomposition (QSD) 
technique [3l| and decomposed the controlled action of 
a two-qubit exp^irH). QSD is known to decompose a 
generic three-qubit gate with the least number of CNOTs 
(20). A minimal number of CNOTs is very important as 
their implementations are orders of magnitude more dif- 
ficult. We found a circuit with 9 CNOTs which is not 
universal in the sense that the decomposition must be 
done for all powers of U individually, or a universal 10- 
CNOT-circuit. The structure of this circuit is shown in 
Figure m The controlled action of nth power of U is sim- 
ply done by multiplication of the angles of Rz rotations 
by n. Details of the decomposition and also all param- 
eters important for a possible experimental realization 
which correspond to the calculations at internuclear dis- 
tance 3.255 ao can be found in the Supplementary In- 
formation. The proposed experiments are undoubtedly 
a challenge for different realizations of quantum compu- 
tation. We regard experimental verification of the usage 
of HF initial guesses in a realistic noisy environment and 
also the performance of both versions of IPEA (A and 
B) proposed in as very interesting. 

Conclusion. - In this work, we have presented the first 
quantum algorithm for 4c relativistic FCI energy compu- 
tations. This algorithm not only achieves an exponen- 
tial speedup over its classical counterpart, but also has 
the same cost (in terms of scaling) as its NR analogue. 
We have proved its functionality by numerical simula- 
tions of calculations of the spin-orbit splitting in SbH. 
We have also proposed and designed the first small-scale 
experimental realizations of relativistic qFCI computa- 
tions. Our algorithm can be used stand-alone or as a 
subroutine of a property algorithm of Kassal et. al. [11 1 
e.g. for calculations of NMR properties. 

This work has been supported by the GACR 
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Figure 5. Scheme of a circuit corresponding to CAS(4,3) calculations on SbH. Empty squares represent generic single-qubit 
gates. _Rz gates are without angle specification. For derivation, details, and all the parameters, see Supplementary Information. 
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SUPPLEMENTARY INFORMATION 



Size of 4c relativistic FCI eigenvalue problem 



Relativistic Hamiltonian approximations 

Our work is based on the 4-component electronic 
Dirac-Coulomb Hamiltonian which in atomic units is 
given as 



In this section, we compare dimensions of non- 
relativistic and 4c relativistic Hamiltonian matrices. In 
the NR case, the Hamiltonian matrix is block diagonal 
according to Ms- Thus for a closed shell system with n 
electrons in m orbitals, the number of determinants is 



^ 1 
H [c{a., ■ Pi) + P'imc^ - (j)^^^] +^ — + Vnn- 

2—1 i<j ^-^ 

(SI) 

We work within the Born-Oppenheimer clamped nu- 
clei approximation which allows to factorize out time- 
dependence of the one-electron problem in the nuclear 
frame. The one-electron operator of the electronic Hamil- 
tonian is accordingly given by the Dirac Hamiltonian in 
the electrostatic potential (/)„„c of clamped nuclei. The 
relativistic energy scale has been aligned with the non- 
relativistic one by subtraction of the electron rest mass. 

The full Lorentz-invariant two-electron interaction can 
not be written in a simple closed form, so approximation 
and thus loss of strict Lorentz invariance is in practice 
unavoidable. In Coulomb gauge the zeroth-order 0{c^) 
operator is given by the Coulomb term employed here. 
This resulting Hamiltonian covers the major part of the 
spin-orbit interaction, including two-electron spin-same 
orbit, as well as scalar relativistic effects. Experience 
shows that the Coulomb term is enough for most chemical 
purposes [H], but for highly accurate molecular spectra 
the Breit (Gaunt) term, carrying spin-other orbit inter- 
action, is recommended. 

A fundamental conceptual problem is that the Dirac- 
Coulomb (-Breit) Hamiltonian has no bound solutions 
due to the one-electron negative-energy continuum solu- 
tions generated by the Dirac Hamiltonian [s^ . We adopt 
the no-pair approximation (NPA), widely used in rela- 
tivistic quantum chemistry [22], in which the N-particle 
basis of Slater determinants is constructed from positive- 
energy bispinors only. This procedure in fact neglects all 
QED effects, but it is justifiable at the energy scale rele- 
vant to chemistry. In particular, the Born-Oppenheimer 
approximation is expected to have larger impact than the 
neglect of QED effects. 

We finally note that the Fock space approach to in- 
clude positronic states within the Dirac-Coulomb(-Breit) 
Hamiltonian approximation [li, 35| should be tractable 
on a quantum computer as well, since the direct map- 
ping (including qubits for positrons) covers the whole 
Fock space generated by a finite basis set. For further 
discussion of the Dirac-Coulomb approximation and how 
to possibly go beyond it the reader may consult Refs.jsJ 
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(S2) 



The relativistic Hamiltonian mixes determinants with 
different Mk values and therefore 



2m 



(S3) 



Using Stirling's approximation in the form 



In to! w —In (27rTO) + toIu m — m for to — ^ oo, (S4) 

and setting m ~ k ■ n, the ratio between the relativistic 
and non-relativistic number of determinants is given by 
the expression 
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Controlled-U circuit design 

In this section, we construct a quantum circuit which 
corresponds to the controlled action of powers of U ~ 
e"^ (see Figure 1 of the paper) for a CI space of di- 
mension 3. For this case, we need two qubits to encode 
the quantum chemical wave function and U has a block 
diagonal structure with 3x3 block of an exponential of 
a Hamiltonian and unity on a diagonal to complete the 
vector space of two qubits. 

We use the Quantum Shannon Decomposition tech- 
nique of Shende et. al. |3l| . It turns out to be very useful 
to generalize the concept of controlled gates to quantum 
multiplexors. A quantum multiplexor is a quantum con- 
ditional which acts on target qubit(s) in a different way, 
according to the state of select qubit(s). If the select 
qubit is the most significant one, then it has the follow- 
ing matrix form 



(S6) 



It performs Uq on the target qubit if the select qubit is 
0) and Ui \i the select qubit is A controlled gate 
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is a special case where Uq = /. More generally, if U is 
a quantum multiplexor with s select qubits and t target 
qubits and the select qubits are most significant, the ma- 
trix of U will be block diagonal, with 2* blocks of size 
2* X 2*. 

A controlled 2-qubit U {c-U2q) is a special case of mul- 
tiplexed U and can be decomposed in the following way 

[ail 



(S7) 



u 




A multiplexed z-rotation in the middle of the circuit on 
the right-hand side (at this stage without angle specifi- 
cation) is in fact a diagonal matrix with second half of 
a diagonal equal to a Hcrmitian conjugate of the first 
one. The circuit equation (|S7[) corresponds to the matrix 
equation 
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W 



(S8) 



Note that right in the equation means left in the circuit 
as the time in a circuit flows from the left to the right. 
We then have 



/ = VDW, 



(S9) 
(SIO) 
(Sll) 



V2 = 
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- 001 + 4>w - 4>ii)- 



be 
has 



extracted 
the form: 



Individual (/)'s in (jS14p can 
from the diagonal of Z?, which 
diag(e-*'^™ ,6-*'^" ,e-^*^« ,e-'''^" ) . 

We would like to emphasize that this is not intended 
to be a decomposition technique for general C/'s, as it 
itself requires classical diagonalization [of U\ see (jSll|) ]. 
A general efficient decomposition of an exponential of a 
Hamiltonian to elementary gates is known only for the di- 
rect mapping [iJl . But this mapping is not suitable for 
small scale experiments due to the relatively high number 
of required qubits and operations thereon. Our aim was 
in fact to prepare the ground for a first non-trivial (more 
than one qubit in the quantum chemical part of the reg- 
ister) experimental realization of (relativistic) quantum 
chemical computation on a quantum computer. 

Because V belongs to the group 0(4) (matrix of eigen- 
vectors of a symmetric matrix) , it can be decomposed 
using only two CNOT gates [sg: 
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H and S are standard Hadamard and phase gates and 
A, B are generic single-qubit gates that can be further 
decomposed e.g. by Z-Y decomposition [l[ 



A single- multiplexed Rz gate (with angle (f>Q for |0) 
state of a select qubit and 0i for |1)) can be implemented 
with the following circuit 



(S12) 
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since ax gates on both sides of Rz turn over the direction 
of the Rz rotation. If we use this approach for demulti- 
plexing the Rz gate in (jS7p , we end up (after some simple 
circuit manipulations) with the following circuit for c-U2q 
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(S13) 
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where 



e'''Rzif3)Ryi-f)RziS). 



(S16) 



There is a highlighted swap gate in (|S15[) which should 
be applied only if the determinant of V is equal to — 1 

The matrix W, on the other hand, is not real as it is 
equal to D^V'^ ([S9|) and can be implemented using three 
CNOT gates (see e.g. [li, SO]). The total count is thus 
9 CNOTs. 

The disadvantage of the aforementioned scheme is that 
W must be decomposed for each power of U individually. 
If we separate W to and D^, is the same for all 
powers of U (eigenvectors don't change) and can be 
up to a non-measurable global phase implemented with 
the following circuit 



-B 



Rzi^e) 



(S17) 



Rz{-^)^Rzm 



Rzifr) 
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Ground state (0^ 


rjxcited state (ij 


900 


-i.UiD4zZ ( o 


-l.UUDOD / DO 


001 


-U.Doo ( 45io 


-U.iooy / yz4 


010 


0.d9d57z37 


-0.39129153 


011 










0.73125768 


-0.00680941 


7 


-0.10311594 


2.21832498 


5 


-0.12107336 


-3.13494247 


A-Eshift 


-6477.89247780 


-6477.89247780 



Table SI. Circuit parameters: rotation angles 0ij, i,j G {0, 1} 
(|S14IS18|) . Z-Y decomposition parameters of A, B (|S15|) and 
energy shifts (core energy + nuclear repulsion) for CAS(4,3) 
calculations of 0'^ and 1 states. For the details see preceding 
text. 

where 



V'S = ^('^OO - 001 - 010 + 01l), 

(yfe = ^(-000 - 001 + 010 + 0ii), (S18) 

^7 = ^(-000 + 0Ol)- 

The circuit for V'' is the same as for V (jSlSp . merely A 
is replaced by and i? by A^. 

Presented 10-CNOT-circuit is universal for all powers 
of U. The only thing one has to do is to multiply the 
angles of Rz rotations in (|S13p and (|S17p according to 
the power of U, e.g. by 2 for the second power. 

Table [511 summarizes the circuit parameters for ground 
as well as excited state calculations described in the pre- 
ceding text. Notice that 0ii is zero in both cases by con- 
struction. To complete the vector space of two qubits, 
we in fact added one eigenvalue of the Hamiltonian equal 
to zero. Other simplification, which originates from the 
block diagonal structure of U, is that A and B matrices 
in the decomposition of V (jSlSP differ only in a global 
phase. Because the global phase is not measurable, we 
present just the angles of rotations. Also only the pa- 
rameters corresponding to A and B are shown. Going to 
their Hermitian conjugates means swapping of /3 and S 
and changing the sign of all of them. 

For the excited state, the determinant of V is equal to 
— 1 and therefore the swap gate in (jSlSp should be ap- 
plied. Because we took Hamiltonian matrices from the 
DIRAC program [i^, the parameters in Table ISll refer to 
the difference between the total energy and core energy + 
nuclear repulsion (A_Bshift). The presented method with 
the parameters form Table [ST] implements the exponen- 
tial e'"^^, as was already mentioned. But in our version 
of the algorithm Q, we in fact need e"*"^^. The obtained 
energy therefore corresponds to the negative of the en- 
ergy. For the negative, the energy guesses ii^max = 3.5 
and -Emin = 2.0 corresponding to the maximum and min- 
imum expected energies were used. 



We don't give any explicit proof that the Quantum 
Shannon decomposition is optimal in the number of 
CNOT gates for the specific case of block diagonal c-U2q- 
However, this conjecture is supported by the fact that we 
also implemented the Group Leaders Optimization Algo- 
rithm (GLOA) of Dashkin and Kais [41| and unsuccess- 
fully tried to find a better circuit (in terms of number of 
controlled operations) with a fidelity error smaller than 
0.01. 



* libor . veis@j h-inst . cas.cz | 

t jiri.pittner@jh-inst.cas.cz, corresponding author 
[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge University Press, 
2000). 

[2] P. W. Shor, Algorithms for quantum computation: Dis- 
crete logarithms and factoring, in Proceedings of 35th 
IEEE Symposium on Foundations of Computer Science, 
pp. 124-134, , 1994, IEEE Press. 

[3] L. K. Groover, Phys. Rev. Lett. 79, 325 (1997). 

[4] R. P. Feynman, Int. J. Theor. Phys. 21, 467 (1982). 

[5] D. A. Lidar and H. Wang, Phys. Rev. E 59, 2429 (1999). 

[6] A. Aspuru-Guzik, A. D. Dutoi, P. J. Love, and M. Head- 
Gordon, Science 309, 1704 (2005). 

[7] H. Wang, S. Kais, A. Aspuru-Guzik, and M. R. Hoff- 
mann, Phys. Chem. Chem. Phys. 10, 5388 (2008). 

[8] J. D. Whitfield, J. Biamonte, and A. Aspuru-Guzik, Mol. 
Phys. 109, 735 (2011). 

[9] L. Veis and J. Pittner, J. Chem. Phys. 133, 194106 

(2010) . 

[10] 1. Kassal, S. P. Jordan, P. J. Love, M. Mohseni, and 
A. Aspuru-Guzik, Proc. Natl. Acad. Sci. 105, 18681 
(2008). 

[11] 1. Kassal and A. Aspuru-Guzik, J. Chem. Phys. 131, 
224102 (2009). 

[12] H. Wang, S. Ashhab, and F. Nori, Phys. Rev. A 79, 
042335 (2009). 

[13] N. J. Ward, 1. Kassal, and A. Aspuru-Guzik, J. Chem. 

Phys. 130, 194105 (2009). 
[14] B. P. Lanyon et ai, Nature Chemistry 2, 106 (2010). 
[15] J. Du et ai, Phys. Rev. Lett. 104, 030502 (2010). 
[16] Z. Li et ai, Nature: Sci. Rep. 1, 88 (2011). 
[17] D. Lu et al, Phys. Rev. Lett. 107, 020501 (2011). 
[18] 1. Kassal, J. D. Whitfield, A. Perdomo-Ortiz, M. H. Yung, 

and A. Aspuru-Guzik, Annu. Rev. Phys. Chem 62, 185 

(2011) . 

[19] D. S. Abrams and S. Lloyd, Phys.Rev.Lett. 83, 5162 
(1999). 

[20] R. B. Grifliths and Chi-Sheng Niu, Phys. Rev. Lett. 76, 
3228 (1996). 

[21] B. A. Hess and C. M. Marian, in Computational Molecu- 
lar Spectroscopy, edited by P. Jensen and P. R. Buenker, 
pp. 169-219, Wiley, Sussex, 2000. 

[22] K. G. Dyall and K. Faegri, Introduction to Relativis- 
tic Quantum Chemistry (Oxford University Press, New 
York, 2007). 

[23] P. Jordan and E. Wigner, Z. Phys. A 47, 631 (1928). 

[24] K. Balasubramanian, Chem. Rev. 89, 1801 (1989). 

[25] T. Fleig, J. Olsen, and L. Visscher, J. Chem. Phys. 119, 



9 



2963 (2003). 

[26] DIRAC, a relativistic ab initio electronic structure pro- 
gram, Release DIRAC08 (2008), written by L. Visscher, 
H. J. Aa. Jensen, and T. Saue, with new contributions 
from R. Bast, S. Dubillard, K. G. Dyall, U. Ekstrom, 
E. Eliav, T. Fleig, A. S. P. Gomes, T. U. Helgaker, 
J. Henriksson, M. Ilias, Ch. R. Jacob, S. Knecht, P. Nor- 
man, J. Olsen, M. Pernpointner, K. Ruud, P. Salek, and 
J. Sikkema (see littp://dirac. chem.sdu.dk). 

[27] L. Visscher and K. G. Dyall, At. Data Nucl. Data Tables 
67, 207 (1997). 

[28] E. Fahri et al, Science 292, 472 (2001). 

[29] I. Buluta and F. Nori, Science 326, 108 (2009). 

[30] R. Gerritsma et al, Nature 463, 68 (2010). 

[31] V. V. Shende, S. S. Bullock, and I. L. Markov, IEEE 
Trans, on Computer- Aided Design 25, 1000 (2006). 

[32] O. Visser, L. Visscher, P. J. C. Aerts, and W. C. Nieuw- 
poort, Theor. Chim. Acta 81, 405 (1992). 

[33] G. E. Brown and D. G. Ravenhall, Proc. Roy. Soc. Lon- 



don A 208, 552 (1951). 
[34] T. Saue and L. Visscher, Four-component electronic 
structure methods for molecules, in Theoretical Chem- 
istry and Physics of Heavy and Superheavy Elements, 
edited by S. Wilson and U. Kaldor, p. 211, Kluwer, Dor- 
drecht, 2003. 

[35] W. Kutzelnigg, Chem. Phys. 395, 16 (2012), DOT: 

10.1016/j.chemphys.2011.06.001. 
[36] T. Saue, ChemPhysChem 12, 3077 (2011). 
[37] W. Liu, PCCP 14, 35 (2012). 

[38] J. Derezihski, lAMP News bulletin January, 11 (2012). 
[39] F. Vatan and C. Wilhams, Phys. Rev. A 69, 032315 
(2004). 

[40] V. V. Shende, L L. Markov, and S. S. Bullock, Phys. 

Rev. A 69, 062321 (2004). 
[41] A. Daskin and S. Kais, J. Chem. Phys. 134, 144112 

(2011). 



